WUB 00-03 
TAUP 2627-2000 



Large Quark Rotations, Neutrino mixing 
and Proton Decay 

Yoav Achiman*f and Carsten Merten^ 

Department of Physics 
University of Wuppertal 
Gaufistr. 20, D-42097 Wuppertal 
Germany 
f School of Physics and Astronomy 
Tel Aviv University 
69978 Tel Aviv, Israel 



March 2000 



Abstract 

Right-handed (RH) rotations do not play a role in the Standard Model, and only the 
differences of the LH mixing angles are involved in V CKM . This leads to the huge 
freedom in the fermionic mass matrices. However, that is no more true in extensions 
of the Standard Model. For example in GUTs large RH rotations of the quarks can 
be related to the observed large neutrino mixing or in particular, all mixing angles are 
relevant for the proton decay. We present a simple realistic non-SUSY 50(10) GUT 
with large RH and LH mixing and study the corresponding nucleon decay rates. 



* e-mail: achiman@theorie.physik.uni-wuppertal.de 
' e-mail: merten@theorie.physik.uni-wuppertal.de 



1 Introduction 



What is the origin of the fermionic masses? This is one of the open questions in the 
Standard Model (SM). The mass matrix entries are arbitrary parameters of the model 
and only the neutrinos are restricted to be massless. 

One can add conjectures for the structure of the fermionic mass matrices to the SM. Many 
different conjectures are known to give the right masses of the charged fermions and V CKM 
(within the experimental errors) and this is clearly an indication that the mass problem 
is far from being solved. 

The main reason for this large freedom is that right-handed (RH) rotations Q are non- 
observable in the SM. Also the observed left-handed (LH) mixing matrix V CKM involves 
only the differences between the mixing angles of the u-like and d-like quarks and the 
individual mixing can be large. Actually, there is already a strong indication from the 
neutrino sector that large rotations are required |2|, [|. 

The considerable freedom in the mass matrices of the SM does not exist in its extensions. 
We know that the SM must be extended for many reasons if not alone to explain the 
fermionic masses and mixing angles. In particular, in its most popular extension, the 
Grand Unified Theories (GUTs) pi], all mixing angles are relevant, as will be explained 
later. In such a framework one cannot change the "weak basis" without changing the 
physics. 

The aim of this paper is to present a simple realistic SO(10) GUT || model with large 
LH as well as RH mixing angles and to use them to study the consequences for nucleon 
decay. Mixing effects were generally neglected in the conventional proton decay models || 
by assuming that the mixing angles are small. 

Large rotations in the quark sector could be a natural reason for the large mixing observed 
in the leptonic sector in terms of neutrino oscillations. In particular, there is a kind of 
duality between the RH mixing of the quarks and the leptonic LH rotations 0] . 
A predictive model for the fermionic masses must involve a family symmetry which dictates 
the texture of the mass matrices and protects it from getting large radiative corrections. 
We used a global U(1) F [||] in the framework of an SO(W) GUT that will add relations 
between the matrix elements. Our aim was to look for the simplest possible realization of 
a realistic model with large mixing angles. We therefore used the famous 50(10) paper of 
Harvey, Reiss and Ramond H and generalized it to asymmetric matrices. Such matrices 
give usually large LH and/or RH mixing angles. 

We did not use non-renormalizable contributions to the mass matrices a la Froggatt and 
Nielsen [flOf] because this method assumes ad hoc physics beyond the GUT and many new 
particles. In addition, the resulting matrix elements are given there in orders of magnitude 
only. This can explain the hierarchy of the masses but not the light see-saw neutrino 
properties. Those are obtained from a product of three matrices and hence predicted up 
to a factor of [0(1)] 3 which may be quite large filjl . 

Most recent models use SUSY GUTs [12] []. However, the available parameter space of 



low-energy SUSY shrinked recently so much that MSSM is on the verge of loosing its 



"naturality" 14]. At the same time solutions without low energy SUSY have emerged 



quite naturally [15] in superstring and M theory, and also the hierarchy problem can be 



1 A general non-hermitian matrix is diagonalized by a bi-unitary transformation. This means that two 
unitary matrices, one from the left and one from the right, are needed. Those matrices are equal only for 
hermitian (or symmetric) matrices 

2 We shall use it also in a forthcoming paper [ psj . 
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solved adding extra dimensions [16|. We think therefore that it is worthwhile to look for a 
non-SUSY GUT which is consistent with all observed experimental facts. The hope is that 
the fine tuning (hierarchy) problem will be solved in the more fundamental theory. Note 
also that conventional GUTs are relatively simple and give much more reliable predictions 
for nucleon decay than SUSY theories. 

Unification of the gauge coupling constants is obtained using an intermediate breaking 
scale [17] Mj ~ 10 11 GeV. This is very useful because M I is also the right mass scale 
for the RH neutrinos needed for the see-saw mechanism as well as for leptogenesis as the 
origin of the baryon asymmetry [18] and the (invisible) axion window [jHj ]. 
The observed neutrino oscillations teach us about the neutrino masses and mixing angles. 
This is the first evidence for physics beyond the SM. Our claim is that observation of other 
phenomena like RH currents, leptoquarks, baryon asymmetry induced by leptogenesis and 
especially nucleon decay can reveal the unknown mixing angles and reduce considerably 
the freedom in the fermionic mass matrices. 

The plan of the paper is as follows. In section 2 we discuss the symmetry breaking that is 
dictated by the requirement of gauge unification and the Higgs representations needed to 
give the correct fermion masses. The mass matrices and our global U(1) F details are given 
in section 3. In section 4 the numerical solutions for the mass matrices are obtained by the 
use of the renormalization group equations (RGEs), and a fit to the observed properties 
of the charged fermions is elaborated. Three solutions are found which are consistent with 
the observed neutrino anomalies [Q, |3], ^0], ^] (except for LSND p2^). All mixing angles 
for those solutions are obtained explicitly and this allows for the calculation of the nucleon 
decay rates in section 5. Section 6 is devoted to the conclusions. 



2 Symmetry breaking of the 50(10) GUT 

We use an 50(10) GUT || which is broken down to the SM via an intermediate scale Mj 
as follows 

SO(W) ^ Gj ^ G SM ^ SU(3) C U(l) em , (1) 



where the intermediate symmetry group is the Pati-Salam one ||23 1 , Gj = G PS = SU (4) c ® 
SU{2) L ®SU(2) R . 

The breaking at M v is done using the Higgs representations $210 (° r ^54 i n models with 
D parity) while for the breaking at M 7 we use the SM singlet of a i-e. the G PS 

representation (10,l,3)i26- For the masses of the Higgs scalars the "extended survival 
hypothesis" |24j is assumed. 

In view of the representation content of the mass terms 

16® 16 = (10 © 126) symm 120 antisymm (2) 

we give the light fermions masses via the VEVs of (1, 2, 2) 10 / 120 and (15, 2, 2) 120 / 126 . 
The RH neutrino masses of order M R ~ M x will be induced via the VEV of the above 
mentioned (10, l,3)i26- 

The Higgs doublet of the SM is therefore a linear combination of the SU{2) L doublets in 
the representations (1,2,2) and (15,2,2). The exact number of representations needed 
for the fermion mass matrices will be given later when we will discuss those matrices. This 
number is needed however to fix the RGEs used to calculate M T , Mu and a [/ (M (7 ) as well 
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as the values of the mass matrices at M 7 As will be explained in the next chapter, we 
used in the RGEs Ni = number of (1,2,2) = 4 and iVi 5 = number of (15,2,2) = 2 as 
well as Ax = number of (10, 3, l)i26 = and An = number of (10, 1, 3)i26 = 1 for the 
G PS symmetry breaking. 
For the matching conditions at M 7 we took 

a-JiM;) = a -\M I ) + -±- (3) 
a- L \M x ) = a-\M r ) (4) 
ot^{M x ) = b - a ~\M I )- 2 -a- l {M I ) + ±- (5) 

while at M v the conditions are 



a-\M v ) = a-J{M v ) + ±- (6) 

= d£(M ) + ± (7) 

= a-^M^ + i- (8) 

In general the matching conditions between gauge couplings belonging to theories with 
symmetry groups Qj and Qk can be written as 

<*7\M t/u ) - ^ S 2 {G j ) = a~ l (M I/v ) - ±- S 2 (Q k ) (9) 

where S2(Gj) is the Dynkin index of the adjoint representation of the group Gj. Details 
can be found in p5| ], the results are given in Table |l|. We checked also the RGEs for the 



Quantity 


Mi 


ai (Mj) 


a 2 (Mj) 


Value 


6.14 • 10 iU GeV 


(46.00)- 1 


(39.86)" 1 


Quantity 


a 3 (Mj) 


a aa (M x ) 


a 2L {Mj) 


Value 


(31.39)" 1 


(55.85)" 1 


(39.86)" 1 


Quantity 




M v 


a ir {Mj) 


Value 


(31.42)" 1 


1.31 • 10 1(i GeV 


(20.08)- 1 



Table 1: Symmetry breaking scales and gauge couplings for N% = 4 and Ni§ = 2 

intermediate symmetry G 7 = G PS 59 D, where D is the discrete L>-parity which requires 
&2L = &2R between M v and M 7 . In this case however M U (G PS <g) D) = 1.16 • 10 15 GeV, a 
value which would lead to a too fast proton decay. 

3 The mass matrices 

The aim of this paper is to present models with large RH and LH mixing angles of the 
quarks that lead to large mixing of the leptons and to study the predictions of those 

3 Note that due to the quark-lepton symmetry of G PS the SO(10) mass relations are valid also at the 
scale Mj. 
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models for the nucleons decay ratios. To generate the mass matrices we use the method 
of Harvey, Reiss and Ramond || as it was realized in [26] for symmetric mass matrices. 



However, in order to have large mixing angles one should rather use asymmetric mass 
matrices. For asymmetric matrices we need to make certain changes in the above scenario 
and add also the antisymmetric Higgs representation $120 • Practically speaking , we 
shall use a global family symmetry U(1) F or Z n that will dictate the neutrino properties 
in terms of the observed masses and mixings of the charged fermions. This symmetry 



will be chosen in such a way that the predictive Fritzsch texture [27] will be realized. 
However, as it is well known the symmetric version of this texture cannot account for the 
large top quark mass J28|| . As we need anyhow asymmetric mass matrices we shall use 
the asymmetric Fritzsch texture which is also known under the name "nearest neighbour 
interaction" model (NNI) f29|. Namely, 

/ A \ 
M = I B C . (10) 
\ D E ) 

In view of the fact that we are actually mainly interested in the predictions for the nucleon 
decay rates which are not sensitive to the details of CP violation we will use for simplicity 
real mass matrices (]. 

The three fermion families and the different Higgs representations in XQ^k^j transform 
under the global U(1) F as follows 

16j -> exp(ia i 6')16 i (11) 
$ fc -> exv(i/3 k 6)<f> k (12) 

The invariance under U (1) F requires therefore that the (5k must obey Oi+Oj = (3k ■ Hence, 
the fermionic part of the mass matrices has the following quantum numbers: 

(ai + a\ a\ + ct2 ai + 03 \ 
ai + a 2 o 2 + «2 02 + 03 • (13) 
Oi + O3 «2 + "3 «3 + «3 / 

To realize the NNI texture ( |10D one sees that only Higgs representations with the charges 
(3 = a\ + 02,02 + 03 and 03 + 03 can couple to the fermions. Also, we still have the 
possibility to couple one Higgs representation to two different combinations i.e. 01 + 02 = 
2o 3 . 

Taking all this into account the Yukawa coupling matrices (at energies ji > Mj) can have 
the structure 

Y£>= \ x 1 I) j : Y;> 6 - I //, i) j : Y^ 






Yi 2 > = x 2 ; Y{ 2 1 = y 2 \ ; Y?) = z 2 \ (14) 



On top of that, we need at least one $126 with a large VEV in (10,1,3) to break the 
SO (10) gauge symmetry at M I . This VEV will generate also the RH neutrinos masses 

4 This can help to solve the strong CP problem while the observed CP violation in the _R"-decay can 
come from a different origin than the CKM matrix. 
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M R ~ Mj. Actually, as is discussed in App. I we will use only one $126 i n addition to two 
$10 and two $120 to generate the asymmetrical mass matrices. The detailed fits required 
four VEVs in the direction (1, 2, 2) and two in that of (15, 2, 2). Those give also the right 
unification as discussed before. 

In terms of the notation and discussion of App. I we obtained the following expressions for 
the mass matrix elements in terms of the VEVs and Yukawa couplings: 

(M d ) 12 = Xl vf + yi uf + Zl vf (15) 

(M d ) 21 = Xivf+y^f-z^ (16) 

(M d ) 23 = x 2 vf+z 2 {vf+uf) (17) 

(M d ) 32 = x 2 vf -z 2 (vf +wf) (18) 

(M d ) 33 = x^+yiw™ 

= (|M"+(|M" < 19 > 

(M e ) 12 = xwf -Sywf+zivf 

= (M (i ) 12 -4y 1 u;( 1) (20) 

(M e ) 21 = x lV f -S yi uf - Zl vf 

= (M^-4^ (21) 

(M e ) 23 = xwf+^vf-ZG,?) 

= (M d ) 2:i -4z 2 Co^ (22) 

(M e ) 32 = x 2 vf-z 2 {vf-3,^) 

= (M d ) 32 + 4z 2 Co d 2) (23) 

(M e ) 33 = Xivf-Sywf 

= (M^-^W 1 ' (24) 



(m u ; 


)12 


= xiv { ^ 


+ vi^u } 




(m«: 


>21 


= 


+ yrf 


- ^ 


(m u ; 


>23 


= x 2 v^ 




+ ^ 2) ) 


(m u ; 


>32 


= x 2 v (2) 




+ ^ 2) ) 


(m„; 


>33 


= xiv^ 


+ 





(25) 
(26) 
(27) 
(28) 



(M< Dir) )i2 = x^-Sy^+z^ 

= (M u ) 12 -4y 1 u;< 1) (30) 

(M™) 21 = x^-Sy^-z^ 

= (M u ) 21 -4 yi u,« (31) 

(M< Di ") 23 = x 2 v^ + z 2 (v^-3u^) 

= (M u ) 23 -4z 2 c^ (32) 

(M™) 32 = x 2 vV-Z2(v™-3u,™) 

= (M u ) 32 + 4z 2 c^ (33) 
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= (M u ) 33 - 4 W 1 ' (34) 

All those elements are obtained from 14 independent Higgs parameters which are products 
of VEVs and Yukawa couplings. They also fix the matrix elements of the RH neutrino 
Majorana mass matrix (in terms of the VEV of (10, 1,3) that breaks G PS — > G SM ), with 
M R ~ Mj being a quasi-free parameter: 

'0 1 

10 j (35) 







yi 






M^ aj) = M R 


2/1 








= yiM R | 









2/1 J 





2/1/2/1/ 

By diagonalizing the charged fermion mass matrices 

uIm u u r = Mi D \ nlM d n R = M^, 

EiM e E fl = Mi D \ utD L = V CKM , (36) 
and the light see-saw neutrino mass matrix 

M 1 ^ ~ -Mf^M^f^Mrf (37) 

we obtained all masses of the charged fermions and those of the three light neutrinos as 
well as the mixing matrices XJ L , R , D L , fl , E ifl and Nj, []. 

Neglecting phases (as explained before) we have 12 masses and 21 mixing parameters. 
All those fix the proton decay rates and other GUT scale effects (like baryon asymmetry 
induced by leptogenesis [|18|). However, in the framework of the SM (with massive neutri- 



nos) the only observable mixing matrices would be V CKM = uJ,D L and U = e|_N^. This 
gives 18 mixing observables []. 

4 Numerical solutions for the mass matrices 

The first step will be to calculate numerically the mass and mixing matrices in (36). This is 



done using the 50(10) relations ([l5| - 04]) for the mass matrices at \x = M 1 = 6.14- 10 10 GeV. 
On the RHS of ([}(]) we use the calculated masses of the charged fermions (see Table (2j) 
and the (real) CKM matrix at fi = M I . Those values are obtained using the RGEs as is 



described in detail in [25]. We then run the mass matrices from M I to M z , diagonalize them 
there and compare the obtained masses and mixings with their experimentally observed 
values. 

Now, without using the neutrino sector, only 13 of our 14 parameters are independent 
because only the combination Z2{vu ) + &u ) appears in the charged fermion equations. 
We have in (^TJ) 30 nonlinear equations. On the other hand each of the (real) 6 mixing 
matrices is parametrized by 3 angles, so altogether there are 18 mixing angles. With the 
13 Higgs parameters we have 31 "unknowns". One of them must therefore be "given" to 
be able to search for solutions numerically. For the given quantity we use the ratio 2/1/2/1 
which determines the RH neutrino mass matrix (^) except for a global factor and we vary 

There is only one neutrino mixing matrix as ~M.^ sht is always symmetric. 
6 Note however, given the mass matrices we have definite predictions for the proton decay, the baryon 
asymmetry and RH currents which can be measured in principle by future experiments. 
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Mass 


m u {Mj) 


m d {Mj) 


m s (Mr) 


Value 


1.16 MeV 


2.38 MeV 


47.4 MeV 


Mass 


m c (Mj) 


771ft (Mj) 


m t {Mj) 


Value 


337.6 MeV 


1360 MeV 


101.2 GeV 


Mass 


m e (Mj) 




m T (Mj) 


Value 


513 keV 


108.14 MeV 


1838.3 MeV 



Table 2: Fermion masses at M T 



Solution 


MSW effect 


yi/yi 


z 2 Cj ( u ] (MeV) 


Model 1 


large mixing 


19 


-3450 


Model 2a 


small mixing 


24 


-9175 


Model 2b 


small mixing 


-18 


-9925 



Table 3: Values of y\jy\ and z^&u for the three representative solutions 



its value between 1 < \y\/yi\ < 1000. By studying the equations in detail one can see 
that, once the value of y\jy\ is given, the neutrino sector of our model is uniquely fixed 
up to the two parameters M R ~ M x and £2^4 • We therefore look if there are solutions of 
the model for reasonable values of yi/yi, M R and Z2&U 1 which predict neutrino properties 
lying in the range allowed by oscillation experiments H, 0, 20, 21]: 



|(U)i 3 | < 0.05 (38) 

0.49 < |(U) 23 | < 0.71 (39) 

0.03 < |(U)i 2 | < 0.05 (small angle MSW) 

or 0.35 < |(U)i 2 | < 0.49 (large angle MSW ) (40) 

50 < AmL/AmL = ( m * - m ^)/(mJ 2 - m 2 ui ) < 1000 (41) 



The value of y\M R will be fixed at the end to give the exact absolute scale of the neutrino 
masses. 

Given these three parameters the model predicts the three neutrino masses and three 
lepton mixing angles in U. We found two regions in the y\jy\-Z2&u ) parameter space 



which obey the atmospheric neutrino requirements together with small angle MSW [30] 
explanation for the solar neutrino puzzle (models 2a, b) and one with the large angle MSW 
(model 1) as can be seen in the Figures ^ and ||. For each region we fixed a representative 
solution (see Table [|) and used it to calculate the corresponding neutrino properties. The 
explicit results for the neutrinos are given in Table |j. 

In this way the explicit LH and RH mixing matrices are also fixed. The corresponding 
mixing angles are given in Table |. They are used to calculate the branching ratios of the 
nucleon decays. 
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Figure 1: Solution 1 (large mixing MSW) in the yi/yi-Z2tOu ) parameter space 
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Figure 2: Solution 2a (small mixing MSW) in the y\/yi-Z20Ju ) parameter space 
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Figure 3: Solution 2b (small mixing MSW) in the y\jy\-Z2^>u parameter space 

5 Calculation of the nucleon decay rates 

The partial decay rate for a given process nucleon — > meson + antilepton is expressed as 
follows: 



(42) 



where M l and A4 r are the hadronic transition matrix elements for the relevant decay 
process. I and r denote the chirality of the corresponding antilepton. A l and A r are the 
relevant coefficients of the effective Lagrangian ( |67| ) given in App.II. A, A L and A R are 
factors which result from the renormalization of the four fermion operators as follows: 



A T 



A 



82 



«i(^h -I 

V a^Mj) ) 

r a^jM^ yi / a 2L (M,) 

Va^Mv)) \a 2L {M v )J \a 2R {M v )J \a 2 (Mj) 
a 3 (M z ) \ I / a 3 (m 6 ) \ h / a 3 (m e ) \ A / a 3 (l GeV) \ f 
(-^j) ' W^zV U a (m 6 )J V a 3 (m c ) J 



Using then ||31| 



a 3 (l GeV) = 0.544 , a 3 (m c ) = 0.412 , a 3 (m b ) 



0.226 



(43) 
(44) 

(45) 
(46) 



one obtains 



1.155 , \A t 



1.071 , \A\ 2 = 23.59. 



(47) 
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Parameter 


T T 1 - TV J" 111 

Value in Model 1 


Value m Model 2a 


Value m Model 2b 


(y 1 M R /M I ) ■ m vi 


-0.0245 eV 


-8.73 • 10" 4 eV 


1.16 • 10~ s eV 


( yi M R /Mj) ■ m U2 


0.0876 eV 


0.355 eV 


-0.467 eV 


{ yi M R /Mj) ■ m V3 


-2.402 eV 


-3.031 eV 


4.365 eV 


u 12 


-0.487 


-0.050 


0.051 


aW) 

t>23 


0.205 


0.506 


0.496 


<M 


0.004 


0.003 


-0.003 


m v .Jm Ul 


-3.57 


-406.3 


-401.2 


m U3 /m y2 
( m 2 vz - ml 2 \ 
\ml 9 - ml t J 


-27.43 
815.4 


-8.54 
71.9 


-9.35 
86.4 



Table 4: Masses and mixing angles of the light neutrinos at M z 



\S\ 2 = (^NuciC^i) ^21 ^s) I ^(^1 ~~ ^2) I ^NuciC^i) ^2> ^3)) i s the probability to find two valence 
quarks of the nucleon at one point in space. We used here the value |5| 2 = 0.012 GeV 3 
given in [||. p j = (1 - - xj) with Xj = "i Mc son/ m Nuci is an ^(6) spin-flavour 

symmetry breaking phase space factor. 

The resulting branching ratios are given in the Tables ^ and [§|. Our model also gives the 
total decay rates. The predicted rates in the different models are given in the Tables 
and |9|. Note however that these predictions have a relatively large uncertainty. This was 



estimated by Langacker []33| to be 

A Wf o = 10 ±0 - 7±1 -°i- yr.,. (.18, 



Taking this into account our predicted representative rates have therefore a chance to be 



observed by SuperKamiokande |3J] and ICARUS [35|. 

Yet the essential predictions of the model are actually the branching ratios. Comparing 
our branching ratios with those of the conventional SO (10) (i.e. without large mixings) 
one sees clearly the suppression of the p, n — > e + X channels relative to the channels 
p, n — > fi + X, u c X . In particular p,n — > /i + 7r, v°K are prominant. Note that in SUSY 
GUTs [12] the dominant decays are into final states involving K mesons []. The special 
properties of our model are clearly reflected in the comparison of ratios as is done in 
Table |To|. 



6 Conclusions 

We presented an 5O(10) GUT with a global U(1) F family symmetry that is consistent 
with all experimental observations. The special thing about this model is that it involves 
large mixing angles for the quarks, in contrast with the conventional expectations. The 
large mixing in the lepton (neutrino) sector results therefore naturally. 
Large rotations have considerable consequences for observables outside the Standard Model. 
In particular calculations of the nucleon decay rates require the knowledge of all mixing 
angles. We found rates which are obviously different from the conventional GUTs. Our 
model is however only one simple example for the effects of large mixing angles. We are 

7 This is also the case for non-SUSY GUTs with maximal RH mixings |36|. 
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Parameter 


Value in model 1 


Value in model 2a 


Value in model 2b 


U L12 


-0.627 


-0.540 


0.275 


o( u ) 

U L23 


-0.055 


-0.056 


-0.075 


0L31 


o.ouo 


u.uou 


n, r\r\r\ 

O.OUO 


0R12 


0.005 


0.006 


-0.012 


U R23 


0.049 


0.051 


0.043 


/)(«) 


n nnn 
U.UUU 


n nnn 
U.UUU 


n nnn 
U.UUU 


aid) 

t>L12 


—0.404 


-0.317 


0.498 


a (d) 
t>L23 


-0.025 


—0.024 


-0.041 


n(d) 
&L31 


— O.Olo 


—0.016 


0.012 


a(d) 

"R12 


0.117 


0.152 


-0.092 


"R23 


0.979 


1.011 


0.679 


n(d) 

VR31 


n nnn 
U.UUU 


n nnn 
U.UUU 


n nnn 
U.UUU 


n le) 
t>L12 


-0.018 


-0.015 


0.015 


Me) 
L 23 


0.733 


-0.060 


-0.089 


a(e) 

V L 31 


n nnn 
U.UUU 


n nn i 
— U.UU1 


n nm 
U.UU1 


U R12 


0.262 


0.314 


-0.301 


n(e) 
U R.23 


-0.063 


0.747 


0.561 


n(e) 
VR31 


0.014 


-0.001 


0.002 


oi") 

Vl2 


-0.487 


-0.050 


0.051 


/)(") 
^23 


0.205 


0.506 


0.496 


oM 
P31 


0.004 


0.003 


-0.003 



Table 5: Mixing angles in the three different models 



studying now a SUSY GUT with all possible phases and its consequences also for other 
GUT observables or RH currents. The hope is that by restricting the values of the mixing 
angles, we will be able to reduce the large freedom in the present models for fermionic 
masses. 
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Rates in % 


Rates in % 


Rates in % 


Rates in % 


of the proton 


(no mixing) 


in model 1 


in model 2a 


in model 2b 


P — > 6^7T^ 


00.0 


91 A 

Z1.4 


9^ 1 
ZO. 1 


97 8 
Z / .0 


p — ► e J\ 




O.l 


9 ft 
Z.D 


4.0 


p — > e + r/ 


l.z 


U.o 


u.y 


1 n 
l.U 


p — ► /i + 7T° 




o.o 


0. f 


0.0 


„ , ,,+ tsO 

p — > /i IS. 


O.o 


ft 
Z.O 


u.y 


1 Q 
l.O 


p -> ^+77 




U.o 


n 9 
U.Z 


n 9 

U.Z 


p — ► e + p° 


0.1 


Q Q 

O.O 


Q O 

O.O 


A 9 

4.z 


p — ► e + o-> 


ICQ 


in o 
1U.0 


1Z. 1 


i a n 
14. U 


p — * e 1\ 




n n 
u.u 


n n 

U.U 


n n 
u.u 


p — ► /i + 




1 Q 
l.O 


n n 
u.y 


U.o 


p — > fI + U> 




4.0 


9 Q 

z.y 


9 a 
z.o 


p — > vflT^' 


on o 

oz.o 


zO.O 


Q £ ft 
OO.O 


ZZ. Z 


p — > I/, il 




9 n 
z.u 


9 n 
z.u 


A 1 
4. 1 








u.o 


n 3 

u.o 


P (i 


n 1 

U.l 


1.0 


n 9 

U.Z 


u.o 


c + 

P —> V e. P 


/I Q 

4.y 


o.y 


K A 
0.4 


A 9 
4.Z 


V ' U c ly 




0.4 


0.3 


0.6 


P P + 




1.4 


0.1 


0.0 


p -► Z/fK*+ 


0.1 


0.0 


0.0 


0.0 


p — > z/? 7T+ 




0.1 


0.1 


0.0 


P — > Z^-fT" 1 " 




0.1 


0.1 


0.0 


P — > 




0.0 


0.0 


0.0 


p -► ^ C K* + 




0.0 


0.0 


0.0 


p -»• e+X° 


56.8 


39.4 


45.1 


51.5 


p -► /j+V° 


5.8 


17.0 


10.6 


11.2 


p -► z/ c V+ 


37.4 


43.7 


44.3 


37.2 



Table 6: Partial decay rates Ti/T of the proton for the case of vanishing fermionic mixing 
in comparison with the rates obtained in our solutions 



Quantity 


Value in model 1 


Value in model 2a 


Value in model 2b 


Tp 


2.54 • 10- 35 yr- 1 
3.94 • 10 34 yr 


2.57 • 10" 35 yr- 1 
3.89 • 10 34 yr 


2.83 • 10- 35 yr- 1 
3.53 • 10 34 yr 



Table 7: Total decay rate and lifetime of the proton for the representative solutions 
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Dppav phannpl 

J—/ C \j OjJ CllcXlllld 


Rafpc: in % 

1V<CXIjCQ 111 /u 


Rafpc: in % 

IV-CXIjCO 111 /U 


Rafpc; in % 

l V^cli u Co ill /(J 


Rafpc: in % 

-HjOjIjCo 111 /u 


ot fnp npiitrnn 

\JL U11C 11C \X 11 wl± 


[nn rniYino'l 

\ 11U 1111-A.lllti, / 


m mnH pi 1 

111 111U LLC1 ±. 


m mon pi 

111 111VJU.C1 LuCXi 


m TTinnpl ^n, 

111 111UU.C1 £i kJ 


n — > e + 7r _ 


CO P 


4U.1 


A f* O 

4b. z 


A P P 

49.9 


n — ► p + ir~ 




ICO 

15.8 


1U.4 


1 p p 
1U.U 


n — ► e + p~ 


y.7 


0.2 


7.1 


1.1 


n — > p + p~ 




O A 

LA 


1.6 


1.5 


n — ► i/ c 7r° 

e 


15.1 


12. U 


lb. 4 


IOC 

lz.5 


n — > zAA 




b.8 


/i p 
4.9 


o o 

8.2 


n ^ u^"q 


0.O 


0.4 


U.6 


U.5 


c 
n — > f^vr 




A O 

4.2 


P o 

0.2 


P 1 

0.1 


n — > v~^K 


1.7 


0.1 


1 P 

1.0 


P P 

0.9 


n — > z/^77 




0.2 


p p 

u.u 


P P 

0.0 


n — ► z/ 7 p° 


1 

z.o 


1.8 


2.5 


1.9 


n — ► iA a; 


7.7 


O.l 


O A 

8.4 


C A 

0.4 


fi > 1/ 1\ 




u.z 


u.z 


n a 

U.4 


C 

n — > v,, p 
pi i 




U.O 


U.U 


0.0 


Q 

n — > v u> 




2.1 


n 1 
U.l 


n i 
0.1 




n n 

u.u 


n i 

u. J- 


n n 

u.u 


o n 

u.u 


n — > h>c-K° 




0.0 


0.0 


0.0 


n -► z/^if 

T 




0.7 


0.4 


0.0 


n — > z/ 7 ^ 




0.0 


0.0 


0.0 


n — »■ i/fp° 




0.0 


0.0 


0.0 


n — ► i/^cj 




0.0 


0.0 


0.0 


n -► v°K*° 




0.0 


0.0 


0.0 


n -► e+X- 


72.6 


46.3 


53.3 


57.6 


n -► 




18.2 


12.0 


11.5 


n -► z/ c X° 


27.4 


35.3 


34.7 


31.0 



Table 8: Partial decay rates Ti/T of the bound neutron for the case of vanishing fermionic 
mixing in comparison with the rates obtained in our solutions 



Quantity 


Value in model 1 


Value in model 2a 


Value in model 2b 


r n 


2.72 • 10- 35 yr- 1 
3.68 • 10 34 yr 


2.80 • 10- 35 yr- 1 
3.57 • 10 34 yr 


3.14 • 10- 35 yr- 1 
3.18 • 10 34 yr 



Table 9: Total decay rate and lifetime of the bound neutron for the representative solutions 
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Ratio 


No mixing 


Model 1 


Model 2a 


Model 2b 


V{p -» e+K") 
T(p -» e+vr ) 


n 


145 


104 

KJ. 1U4 


169 

KJ. lUii 


r(p -» / U+7T U ) 

r(p -> ^+ir ) 


o 


3.27 


6.33 


3.11 


r(p -> v c k+) 


0.003 


0.098 


0.064 


0.157 


T{p e + vr u ) 

r(p -» Z/ C 7T+) 


1.040 


0.618 


0.693 


0.993 


T(n — > /-i + vr _ ) 
r(n — > e + 7r~) 


n 
u 




u.zzo 


n 9nn 


T(n -► /x 
T(n — ► e + p~) 





0.387 


0.225 


0.195 


r(n -» ^ c i^ u ) 
r(n — ► ^ c 7r°) 


0.113 


0.469 


0.347 


0.722 


r(n — > e + 7r _ ) 
r(n — > i/ c 7T°) 


4.16 


2.48 


2.78 


3.96 



Table 10: Ratios of some partial nucleon decay rates 

Appendix I: Fermionic masses in 50(10) GUTs 

Dirac masses for the charged fermions are of the form m{^f : ) T L and therefore transform 
under SO (10) as follows: 



16 16 = (10 126) symm © 120 a ntisyn 



(49) 



The S'O(IO) Higgs representations which can give masses to the fermions have the repre- 
sentation content 

10 — (1,2,2)0(6,1,1) (50) 
120 — > (1,2, 2) 0(15, 2, 2) 0(6, 3, 1)©(6,1, 3) ©(10, 1,1) ©(TO, 1,1) (51) 
126 — ► (15,2,2) © (10,3,1) © (TO, 1,3) © (6,1,1) (52) 

under G PS = S , C/(4) C SU(2) L SU(2) R , while for the fermions in the 16 

16 — > (4,2,1)0(4,1,2) (53) 

Hence the Dirac mass terms transform as follows 

(4, 2, 1) (4, 1, 2) = (15, 2, 2) © (1, 2, 2) (54) 

The different Yukawa couplings and the corresponding VEVs of $10, and ^126 are 



given in Table 11 



Taking now the Clebsch-Gordan coefficients into account the mass matrices get the fol- 
lowing contributions 



M e 

M (Dir) 



v d Y 10 + u d Y 126 + (v d + Lb d ) Y 120 
v d Y 10 - 3w d Y 126 + (v d - 3u d )Y 120 

26 20 
i^uY 10 3lo u Y 126 + (t} u 3lD u )Y 12 o 



(55) 
(56) 
(57) 
(58) 
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Higgs representation 


(1,2,2)SJ 


(1,2,2)S 


(15,2,2)^ 


(15,2,2)W 6 


Yukawa coupling matrix 


v (i) 

1 10 


1 120 


1 120 


1 126 


Vacuum expectation values 






~(fc) 





Table 11: Higgs couplings and vacuum expectation values in S0( 10) GUTs 



The RH neutrinos can acquire also Majorana masses in term of the (10, 1, 3) component 

of $126: 



(4, 1,2) (»(4,1,2)® (10, 1,3) 



M 



(Maj) 
UR 



MrY, 



(1,1,1) 



M;Y 12 



(59) 
(60) 



The corresponding VEV Mj is responsible for the symmetry breaking step G PS — ► G SM - 
The neutrinos will have therefore the 6x6 mass martix 



M 







Mr 



(Dir) 



'ML Dir) ) T M^ aj) 



(61) 



Using the fact that the non-vanishing entries of the Majorana mass matrix are much larger 
than those of the Dirac matrix one can approximately block diagonalize the 6x6 matrix 
and obtain the see-saw matrix [0^1 



(62) 
(63) 



M^ ht « -MS Dir) (M^ j) ) -1 (MS Dir) ) T 

as well as 

M heavy _ 

M)) ght is symmetric and therefore can be diagonalized using one unitary matrix N/y 

N^M^IN^ = M" ght(D) (64) 
Neutrino oscillations are induced via the leptonic analogue to the CKM matrix 



U = E[N, 



(65) 



Appendix II: Effective 50(10) Lagrangian for nucleon decays 

The baryon number violating part of the 5O(10) Lagrangian (without fermionic mixing) 
is known to be 01 



C 



= ^=X; (e aPl u c SY< + d La re + L +d Ra re + R ) 



eh 



u 



+ ^= K a ( - ^dVrdl - u La Yv 

+ ^= XI {d La Ye~ L + d Ra re~ + u La YK L + ^h^Kh) 
+ h.c. 



RaH ^r) 



(66) 
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Taking all possible fermion mixings into account one obtains [4C] the effective four fermion 
Lagrangian 



A, 



+ A 3 {e aM u c L -Y<) [Ptl.dl) + A A (e aM u c L -Y<) (UtlA) 

+ A 5 {e aM u c L ^Y<){e + L l,st) + A, (e Q/3 XVX) (^7m*S) 

+ A 7 {e^u^YO (fiiwl) + A 8 (e aM u c SY<) [fl^s%) 

+ A 9 (e a ^rdt){v? R 7,d" R ) + A w (s a ^rdl)(^ R l,d%) 

+ A n {e a ^rdi){v? R ^s a R ) + A 12 {e a ^Ydi){u c ; R ^s R ) 

+ A 13 {e a ^rsl){^ R 7,d R ) + A 4 (^^ 7 7 m ^)(^7m<) 

+ A 15 {e a ^rdl){^ R ^d a R ) + A w (e a01 u c L -<rd f t){^ R 7,s a R ) 

+ A 
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+ ( terms with two s quarks ) 
+ ( terms with c ,6 and t quarks 
+ ( terms with f£ R and v : 



e,(i,TL 



+ h.c. 



where the coefficients A t are given as follows 25]: 
At 



-- G ((U«) n (UOn + (U K ) 21 (U0 2 i + (U«) 3 i(U t ) 31 
•((E B ) n (D L ) n + (E«) 21 (D £ ) 21 + (E„) 31 (D i ) 31 
+ G ((U fl ) 11 (D,) 11 + (U fl ) 21 (D L ) 21 + (XJ R ) 31 (D L ) 31 
■((E R ) U (U L ) U + (E fl ) 21 (U L ) 21 + (E H ) 31 (U,) 31 
A 2 = ^((UHjnCU^n + CUHkCUOai + CU^aiCU^a! 

•((E i ) 11 (D fl ) 11 + (E t ) 21 (D a ) 21 + (E t ) 31 (D a ) 31 
+ G' ((D^nCUiJn + (D H ) 21 (U,) 21 + (D R ) 31 (U L ) 3l 
•((E £ )n(U fl ) n + (E £ ) 21 (U H ) 21 + (E £ ) 31 (U„) 31 
A, = ^((UHjuCU^n + CUH^CUOai + CU^CU^ 
•((E«) 12 (D £ ) 11 + (E fl ) 22 (D L ) 21 + (E H ) 32 (D i ) 31 
+ G ((U fl ) 11 (D i ) 11 + (U H ) 21 (D £ ) 21 + (V R ) 31 (T> L ) 31 
•((E B ) 12 (U L ) n + (E B ) 22 (U L ) 21 + (E a ) 32 CU £ ) 31 
A 4 = G((U«) 11 (U i )u + (U ii ) 21 (U0 2 i + (U Ji ) 31 (U £ ) 31 
•((EJ^D^ + (E L ) 22 (D fl ) 21 + (E L ) 32 (D JJ ) 31 
+ G' ((D^nCU^u + (D fl ) 21 (U,) 21 + (D fl ) 31 (U i ) 31 
•((EOia^n + (E L ) 22 (U«) 21 + (E L )3 2 (U fl ) 31 
A 5 = G((U B ) 11 (U i ) 11 + (U fl ) 21 (U i ) 21 + (U i? )3 1 (U L ) 31 
•((E H ) U (D £ ) 12 + (E fl ) 21 (D L ) 22 + (E fl ) 31 (D L ) 32 
+ G ((U B ) n (D L ) 12 + (U H ) 21 (D,) 22 + (U i? ) 31 (D i ) 32 
■((E«)h(Ul)h + (E B ) 21 (U L ) 21 + (E fl ) 31 (U L ) 31 



(67) 
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G ((U R ) U (U L ) U + (U,) 21 (U,) 21 + (U K ) 31 (U L ) 3 i) 
•((E L ) n (D«) 12 + (E L ) 21 (D fl ) 22 + (E i ) 31 (D«) 32 ) 
& ((D fl ) 12 (U L ) n + (D fl ) 22 (U L ) 21 + (D fl ) 32 (U L ) 31 ) 
•((E L ) n (U K ) n + (E L ) 21 (U«) 21 + (E i ) 31 (U«) 31 ) 
G ((U fl ) u (U L ) u + (U fl ) 21 (U L ) 21 + (U fl ) 31 (U L ) 31 ) 
•((E fl ) 12 (D L ) 12 + (E fl ) 22 (D L ) 22 + (E fl ) 32 (D L ) 32 ) 
G ((U fl ) u (D L ) 12 + (U fl ) 21 (D i ) 22 + (U,) 31 (D L ) 32 ) 
•((E fl ) 12 (U L ) n + (E R ) 22 (XJ L ) 21 + (E R ) 32 (\J L ) 3l ) 
G ((U fl ) u (U L ) u + (U,) 21 (U L ) 21 + (U fl ) 31 (U,) 31 ) 
■((E i ) 12 (D fl ) 12 + (E L ) 22 (D R ) 22 + (E i ) 32 (D fl ) 32 ) 
& ((D J j) 12 (U i ) 11 + (D fl ) 22 (U,) 21 + (D^tUOai) 
■{(E L ) 12 (V R ) U + (E L ) 22 (U J? ) 21 + (E L ) 32 (U R ) 31 ) 
-G ((U S ) 11 (D I ) 11 + (U«) 21 (D L ) 21 + (U fl ) 31 (D i ) 31 ) 
•((N,) 11 (D JJ ) 11 + (N L ) 21 (D«) 21 + (N i ) 31 (D fl ) 31 ) 
-& ((D fl ) n (D,) n + (D fl ) 21 (D L ) 21 + (D,) 31 (D L ) 31 ) 
■((Nl)ii(U s ) u + (N L ) 21 (U«) 21 + (N L ) 31 (U R ) 3l ) 
-G ((U fl ) 11 (D i ) 11 + (U«) 21 (D i ) 21 + (U fl ) 31 (D,) 31 ) 
■((N £ ) 12 (D fl ) n + (N L ) 22 (D fi ) 21 + (N i ) 32 (D ii ) 31 ) 
-& ((T> R ) U (T> L ) U + (D fl ) 21 (D,) 21 + (D K ) 31 (D i ) 31 ) 
■{(N L ) 12 (V R ) U + (N i ) 22 (U J? ) 21 + (N L ) 32 (U fl ) 31 ) 
-G ((U fl ) 11 (D,) 11 + (U fl ) 21 (D L ) 21 + (U fl ) 31 (D i ) 31 ) 
•((Nl)h(D J j) 12 + (N.^D^ + (N i ) 31 (D ii ) 32 ) 
-G' ((D K ) 12 (D L ) n + (D fl ) 22 (D i ) 21 + (D fl ) 32 (D i ) 31 ) 
■((Nl)h(U h ) u + (N L ) 21 (U fl ) 21 + (N i ) 31 (U J j) 31 ) 
-G ((U^nlDJn + (U«) 21 (D,) 21 + (U fl ) 31 (D £ ) 31 ) 
■((N i ) 12 (D fl ) 12 + (N L ) 22 (D fi ) 22 + (N i ) 32 (D fl ) 32 ) 
-& ((D fl ) 12 (D i ) 11 + (D R ) 22 (D L ) 2 i + (D fl ) 32 (D i ) 31 ) 
'((Ni)i 2 (U fl ) 11 + (N i ) 22 (U H ) 21 + (N i ) 32 (U J{ ) 31 ) 
-G {(V R ) U (B L ) 12 + (U,) 21 (D L ) 22 + (U fl ) 31 (D L ) 32 ) 
■((NJnCD^n + (N L ) 21 (D fi ) 21 + (N i ) 31 (D H ) 31 ) 
-G' ((D fl ) 11 (D i ) 12 + (D R ) 21 (D,) 22 + (D B ) 31 (D L ) 32 ) 
■((NJnlUJu + (N L ) 21 (U«) 21 + (N i ) 31 (U J j) 31 ) 
-G ((U fl ) u (D L ) 12 + (U«) 21 (D i ) 22 + (U fl ) 31 (D i ) 32 ) 
■((N i ) 12 (D fl ) 11 + (N L ) 22 (T> R ) 21 + (N L ) 32 (T> R ) 31 ) 
-& ((D fl ) 11 (D i ) 12 + (D R ) 21 (D L ) 22 + (D ii ) 31 (D i ) 32 ) 
■((Nz,)i 2 (U J{ ) 11 + (N i ) 22 (U H ) 21 + (N i ) 32 (U J{ ) 31 ) 
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.4 



15 



A 



16 



A 
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-G ((U,) n (D L ) n + (U fl ) 21 (D L ) 21 + (U H ) 31 (D i ) 31 ) 
•((N,)i 3 (D«) n + (N L ) 23 (D fl ) 21 + (N L ) 33 (D fl ) 31 ) 
-& ((Dj?)n(D i ) 11 + (D fl ) 21 (D L ) 21 + (D fl ) 31 (D L ) 31 ) 
•((N L ) 13 (U fl ) n + (N L ) 23 (U fl ) 21 + (N L ) 33 (U fl ) 31 ) 
-G ((U fl ) n (D L ) n + (U«) 21 (D i ) 21 + (U fl ) 31 (D i ) 31 ) 
•((N L ) 13 (D fl ) 12 + (N L ) 23 (D fl ) 22 + (N L ) 33 (D fl ) 32 ) 
-& ((D fl ) 12 (D L ) n + (D fl ) 22 (D L ) 21 + (D fl ) 32 (D L ) 31 ) 
•((NzOiaCUjOu + (N L ) 23 (U fl ) 21 + (N £ ) 33 (U fl ) 31 ) 
-G ((U fl ) n (D L ) 12 + (U fl ) 21 (D L ) 22 + (U fl ) 31 (D L ) 32 ) 
•((N0i 3 (D«) u + (N L ) 23 (D fl ) 21 + (N L ) 33 (D fl ) 31 ) 
-& ((D fl ) n (D L ) 12 + (D fl ) 21 (D L ) 22 + (D fl ) 31 (D i ) 32 ) 
•((N L ) 13 (U fl ) n + (N L ) 23 (U fl ) 21 + (N £ ) 33 (U«) 31 ) 



2 

x'y'i 



where Mj y 



We used here the definitions G = gf,/2M^ Y and G' = gi/2M^ 
My is assumed. 

The coefficients A { are connected to the hadronic transition amplitudes of the elementary 
processes responsible for the nucleon decays. The independent amplitudes are given in 
tables ^ and 13. 



Decay process 


Lagrang 


ian term 


Amplitude (-\/30) 


Coefficient 


p| -> e + R u c i u[ 
Pi -> e^u c i ui 
pi -> e + R d c i di 
pi -> e + R d c i di 


{uZYuA fe^7 M dJ 
K7"^)(efl7^ s ) 

(«i7"«i.)(ei7/irfH) 
(^Vw J (e^dj 


-4 
-8 
-8 
+2 


A 2 
A 2 
A 2 
A 2 


Pi -> 4w c | «i 
Pi -> e+d c i d| 


fe7^ i )(e+7 M d fl ) 


-10 
-2 


A 2 
A 2 


p| -> e + R s c i di 
Pi -> e+s c i d| 




(e£7M s fl) 
(efl7M s fl) 


-8 
+2 


A 6 
A> 


Pi -» 4* c i d j 


{u c h Yu Ll 




-2 


\ 


pt -» ^ C T «i 
pT -» fSt<*°i «T 


(«i7 M <y (^ e fl7^fl) 

(u^dj (Pf^M^) 


+4 
-10 


A 9 
A 9 


pi ^ik^i wj. 


(^7 M ^) (^7^) 


-8 


A 9 


pT -» ^ C T «i 

Pt^ ^a°i«T 

PT ^eR^i U i 

pi -> «&« c i «T 


(«?7 M sj Kh7,A) 
(u£ 7 'X) Kh7,A) 


+4 
+2 

-12 


Ai 
Ai 
A3 

A\3 


Pi ^efl^i u| 

pi -> f£l« i <4 


(«i7"d i ) (^7m«k) 
(u^7^ L ) (P^dj 


+4 
-12 


A n 

A\3 



Table 12: Decay amplitudes for the elementary processes of proton decays 
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Decay process 


Lagran^ 


dan term 


Amplitude (V30) 


Coefficient 


nf - 


-> e> c T di 


{u c L Yu L 


) (e^T^dfl) 


-4 


^2 


n| - 






) (e^A) 


+10 


A 2 


n{ - 


- e> c I di 


(u c L Yu L 


) {efr/M 


+8 


A 2 


n] — 


> ^> c T«i 


{u c L rd L ) 


[v?n%d R ) 


+8 


A 9 


n] — 


> ^ c fl n c | lit 


{u° L Yd L ) 




-2 


A 9 


n] — 


♦ ^d c T di 


Hrd L ) 




+4 




n] — 


* ^ R d c i d] 


{u c L yd L ) 


(Ve R l^d R ) 


+8 


A 9 


n[ — 


> i/f^W 


mrdA 


(v? R 1vd R ) 


+2 


A 9 


ni — 


» ^d°I dl 


(u c L Yd L ) 


{v° R l»d R ) 


+10 


A 9 


n| — 


♦ ^s c T di 


{u c L rd L ) 




+4 


A n 


n| — 


♦ ^ R s C i d| 


(u c L rd L ) 


{Ve R l»S R ) 


-4 


A n 


n| — 


♦ "f** c t di 




{vf R l»d R ) 





A%3 


n| — 


* ^«s C i d| 




(v? R l»d R ) 
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Table 13: Decay amplitudes for the elementary processes of neutron decays 
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